Elastic and magnetic effects on the infrared phonon spectra of MnF2 
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We measured the temperature dependent infrared reflectivity spectra of MnF2 between 4 K and 
600 K. We show that the phonon spectrum undergoes a clear renormalization at Tn- The ab-initio 
calculation we performed on this compound accurately predicts the magnitude and the direction of 
the changes in the phonon parameters across the antiferromagnetic transition, showing that they 
are mainly induced by the magnetic order. In this material, we found that the dielectric constant 
is mostly from phonon origin. The large change in the lattice parameters with temperature seen 
by X-ray diffraction as well as the A2n phonon softening below Tn indicate that magnetic order 
induced distortions in MnF2 are compatible with the ferroelectric instabilities observed in Ti02, 
FeF2 and other rutile-type fluorides. This study also shows the anomalous temperature evolution 
of the lower energy ±L/u mode in the paramagnetic phase, which can be compared to that of the Big 
phonon seen by Raman spectroscopy in many isostructural materials. This was interpreted as being 
a precursor of a phase transition from rutile to CaCb structure which was observed under pressure 
in ZnF2. 

PACS numbers: 63.20.kk, 78.30.-j, 63.20.dk, 75.85.+t, 63.20.-e 



I. INTRODUCTION 



In magnetoelectric multiferroic materials, ferroelectric- 
ity coexists with a magnetic order. Based on the origin 
and strength of the coupling between the ferroelectric and 
the magnetic order parameters, these materials can be di- 
vided up into two general classes.^ In type-I multiferroics, 
such as BiFeOsf ferroelectric and (anti) ferromagnetic 
transitions are independent and weakly coupled. In 
type-II multiferroics, such as TbMnOs^ ferroelectricity 
is a consequence of the magnetic ordering and a strong 
magneto-electric coupling is present. In the latter com- 
pounds, the interaction between magnetic ordering and 
the lattice may generate the structural distortions leading 
to the appearance of a permanent electrical polarization. 
Although the microscopic origins for this magnetoelec- 
tric coupling are still under debate, the appearance of 
an electrical dipole moment has important consequences 
on the polar, infrared active, phonon spectra. Very few 
studies of the phonon changes in these type-II multifer- 
roics exist to date. Nevertheless, Schmidt et al.^ showed 
that in TbMnOa phonon parameters are renormalized by 
about 1-2% and that the changes are linked to T/v rather 
than the ferroelectric transition. 

Before tackling the structurally complex type-II mul- 
tiferroics, it is interesting to see the effect of magnetic 



ordering on the infrared phonon response. Several in- 
frared studies on spinels, showed that the magnetic 
ordering leads to infrared phonon splitting. However, 
the phonon splitting has systematically been attributed 
to the magnetic frustration present in these compounds, 
a vision supported by ab-initio calculations Previous 
studies on antiferromagnetic MnO (Ref. TO) and CoO 
(Ref. 11) showed a strong phonon renormalization at T/v- 
However, both materials also undergo a structural phase 
transition from cubic to rhombohedral (MnO) or tetrago- 
nal (CoO) at Tn. These phase transitions make it harder 
to separate what is the pure magnetic effect on phonons 
from what is due to the structural symmetry change. 

In this perspective, manganese fluoride (MnF2) 
emerges as a system of choice, once it is a very 
well characterized commensurate antiferromagnet with 
Tn — 68 K.^^ Its simple paramagnetic rutile structure^^- 
(P42/mnm or D^^) remains the same below Tat, where 
the spins align antiferromagnetically along the 1^4 axis. 
Hence effects on the phonon spectra have a magnetic ori- 
gin. 

Temperature dependent Raman spectra of MnF2 show 
phonon frequency changes at Tat.^^ As this material has 
an inversion center, Raman active phonons are not in- 
frared active and vice versa. Therefore, detailed knowl- 
edge of infrared phonons in MnF2 is paramount to grasp 
changes that would affect dielectric ordering below Tn- 
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We expect that these data on MnF2 help to set a base- 
Hne to understand the effect on phonons of the magnetic 
transitions in type-II multiferroic materials. 

To date, only the room temperature infrared phonon 
spectra have been measured for MnF2.^^ According to 
dilatometric measurements^^ MnF2 lattice parameters 
are renormalized across the antiferromagnetic transition. 
This magnetostrictive effect is the hallmark of the cou- 
pling between dipole excitations and magnetic order in 
this compound. 

Although MnF2 is a classical antiferromagnet, it is 
closely related to multiferroic materials. Its isostruc- 
tural compound Ti02 is a quantum paraelectric (in- 
cipient ferroelectric) as determined by infrared, Ra- 
man and dielectric measurements,^^ as well as ab-initio 
calculations.^^ Similar responses, in particular renormal- 
ization of the phonon spectra, were observed in other 
rutile fluorides such as FeF2,^^ ZnF2,^^ and NiF2.^^ 

In this paper we show a detailed, temperature- 
dependent infrared study of the phonon spectra in MnF2 
along the a6-plane and the c-axis. Our results show that 
the phonon spectra have marked changes at Tat. The in- 
frared data are complemented by low temperature X-ray 
diffraction and ab-initio calculations. First principles cal- 
culations predict the proper phonon frequencies for both 
directions and find the correct frequency shifts at T^. 
Our results show that the dielectric constant of MnF2 is 
mostly from phonon origin. The large change in the lat- 
tice parameters with temperature and phonon softening 
in the antiferromagnetic phase suggest that MnF2 dis- 
tortions are compatible with the ferroelectric instabilities 
observed in Ti02,^^ FeF2,^° and ZnF2^ 



II. METHODS 

The MnF2 single crystal used in this experiment was 
grown by the Czochralski method. X-ray and infrared 
measurements were done on different platelets cut from 
the same bulk and containing ac and ah planes. The 
ac sample for optical measurements was polished with 
a 15° wedge to avoid interference fringes from the back 
surface reflectance. The measured faces for both opti- 
cal samples were polished with Ifim diamond powder to 
a mirror like surface. Typical sample surface sizes are 
5x5 mm^. Their antiferromagnetic transition was mea- 
sured on a Quantum Design MPMS-5 squid magnetome- 
ter with a magnetic field of 1000 gauss. Figure [T] shows 
the magnetic susceptibility (x) measured with the field 
applied parallel and perpendicular to the c axis. The 
H \\ c curve shows a classical antiferromagnetic ordering 
behavior where x decreases below Tn reaching a vanish- 
ingly small value at T = K. This indicates a negligible 
amount of impurities in this sample. When ± c, x is 
dominated by a spin canting response below T^. 

We measured the infrared reflectivity spectra near nor- 
mal incidence at 31 different temperatures between 4 K 
and 300 K, on the ac plane sample with the electric field 
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FIG. 1: (color online) Magnetic susceptibility of MnF2 mea- 
sured with H = 1000 gauss. The solid (blue) line is the data 
for H J_ c and the dashed (red) line is obtained for H \\ c. 
The diagram represents the lattice and magnetic structure 
of the rutile tetragonal phase, where a — b — 4.874 A and 
c = 3.300 A. 



of light parallel and perpendicular to the c direction (D4 
axis). To determine the absolute reflectivity, we used 
an in situ gold overfilling technique^^^ The sample is at- 
tached to a cold finger of an ARS Helitran cryostat and 
we measure its reflectivity with respect to a reference 
stainless mirror for all temperatures. We then evaporate 
a thin layer of gold on the sample producing a mirror 
with the same surface quality and shape as the sample. 
The reflectivity of this gold mirror is measured, at sev- 
eral temperatures, against the same stainless steel refer- 
ence and the infrared spectrum is obtained by dividing 
the uncoated sample by its coated version. Finally, we 
multiply the result by the absolute reflectivity of gold 
to get the absolute reflectivity of the sample. The ac- 
curacy of the absolute reflectivity is better than 1% and 
the relative error between different temperatures is of the 
order of 0.1%. The far infrared (10-700 cm~^) data were 
collected with a Bruker IFS113v interferometer. Higher 
frequency spectra (500-7500 cm~^) were obtained with a 
Bruker IFS66v spectrometer. In the overlapping region 
the spectra agree within 0.5%. 

We also measured the reflectivity of the ab sample, in 
the (100-700 cm~^) range, at high temperatures (300- 
600 K) with a TS-1500 Linkam hot stage. No polarization 
was used in these measurements as a and b are equivalent 
axes. An Al mirror served as a reference. As no in situ 
gold evaporation was used above room temperature, the 
data were corrected so that 300 K measurements in the 
cryostat and in the hot stage match. 

X-ray diffraction measurements were carried out on 
the National Synchrotron Light Source beamline X21. 
A Si(lll) double-crystal monochromator was used to set 
the incident energy of the beam to 11.5 or 13.5 keV, and 
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a Pt-coated mirror focused the beam down to ~ 1 (ver- 
tical) by 2 (horizontal) mm^. The MnF2 crystals were 
inserted in an Oxford superconducting magnet, which is 
mounted on a 2-circle diffractometer with a horizontal 
scattering geometry. A LiF(200) analyzer was used, and 
two parallel reflections were measured for both ac and ab 
plane samples to obtain the a(= b) and c lattice param- 
eters. 

These measurements were complemented with ab- 
initio density functional (DFT) calculations. In order to 
study the effects of the magnetic order on the phonon 
parameters, we performed the geometry optimization 
and phonon calculations for two different spin config- 
urations. The first one is a true ferromagnetic (FM) 
configuration with all spins aligned. The second con- 
figuration is a pseudo-antiferromagnetic (AFM) configu- 
ration where nearest neighbor spins are anti-aligned, as 
shown in Fig. [1] (up-down configuration). This configu- 
ration is only pseudo-antiferromagnetic since it is not an 
eigenfunction of the total spin operator. The true an- 
tiferromagnetic state is a superposition of the up-down 
and down- up configurations, as well as quantum fluctua- 
tions on them. Such a correct calculation is unfortunately 
not feasible for an infinite system. The calculations were 
done using the CRYSTAL06 package*^ Three different 
functional were used for the calculations, namely LDA, 
and two hybrid functional BSLYP^^ and BIPW.^^ An 
atomic basis set of valence 2-( quality and small core 
pseudopotentials^^ were used for the Mn^+ ions and an 
all electrons basis of 3-^ quality was used for the F~ 
ionsj-^^ As usual, the LDA functional underestimates the 
lattice parameters by a few percent while the two other 
functionals slightly overestimate them. The best fit is 
reached for the BIPW functional with an error of 0.2% 
in the a, 6 direction and 1% in the c direction. All ab- 
initio results further presented in the paper will thus refer 
to the BIPW functional. 
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III. RESULTS 

The solid lines in the top panel of Fig. [2] show the in- 
frared reflectivity at selected temperatures with the elec- 
tric field of light lying on the a6-plane. These spectra 
show three clear phonon peaks. The solid lines in the bot- 
tom panel are the reflectivity for ^ || c and show a single 
phonon along this direction. The infrared modes found 
agree with group theory predictions. Indeed, the irre- 
ducible representation decomposition for the rutile struc- 
ture is Aig © A2g © 2A2u © Big © 2Biu © B2g © © 4.Eu. 
Four of these modes are Raman active — Aig^ Big^ B2g 
and Eg — and four are IR active — 3Eu © A2u • The Eu 
modes are xy degenerate and represent the ab plane spec- 
tra. From low to high frequencies, we define these modes 
as Eui, Eu2, and Eus. The A2u mode has z symmetry 
and is the only phonon along the c axis. 

To extract quantitative information from these data, 
we simulated the spectra with a multiple Lorentz oscilla- 



FIG. 2: (color online) Reflectivity with light polarized in the 
a5-plane (top panel) and along the c-axis (bottom panel) for 
MnF2. The solid lines are experimental results and the dot- 
ted lines are Lorentz fits using Eq. (TJ In each panel, the 
dashed line is the reflectivity predicted by ab-initio calcula- 
tions. These curves were generated by plugging into Eq.[T]the 
first principles values for Qto and Ae and using the 7 val- 
ues obtained from the experimental fits at 5K. The ab-initio 
values come from the AFM calculation. 

tor model for the dielectric function: 

where £00 is the contribution from electronic transitions 
to the dielectric function, and each phonon is described 
by a resonance frequency ^rOk^ oscillator strength 
Ae/e, and damping 7/^. The reflectivity at normal inci- 
dence is given hy R = \1 — ^/e\'^/\l + V^p. 
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Typical fitting results are shown as dotted lines in both 
panels of Fig. [2l Overall, the model reproduces the data 
very well. However, modes Eus and A2u^ as well as Eui 
at high temperatures, have structures that cannot be re- 
produced by the Lorentz model. This structure is likely 
related to a breakdown of the harmonic approximation^^ 
and/or two-phonon absorption, even though we cannot 
rule out the presence of a small symmetry breaking lat- 
tice distortion. Nevertheless, in general, the parameters 
obtained through a Lorentz modeling of the data assum- 
ing a D\f^ symmetry are a very good first approximation. 

To ascertain that this is the case with our data, we also 
used Kramers-Kronig transformations, which are model 
independent. These transformations require knowledge 
of the reflectivity in the full spectral range, from zero 
to infinity. As our measurement range is limited, we 
extrapolated the low frequency range as a constant re- 
flectivity. For high frequencies we used a constant up to 
80 000 cm~^ followed by a free-electron approximation 
{R (X uj~^). The 80 000 cm~^ limit was chosen to avoid 
unphysical negative values (albeit within error bars) in 
the imaginary part of the dielectric function {£2)- 

In Fig. [3] we compare Kramers-Kronig results to 
Lorentz fit parameters. The left panel shows £2 for the 
Eui mode. For clarity, each curve was normalized by 
its maximum value and shifted vertically by an amount 
proportional to its temperature. The peak in £2 hap- 
pens at the phonon resonance frequency. The squares 
are the flro frequencies obtained from the Lorentz fits 
to the data. Although a small difference in frequency is 
seen, both methods give the same temperature depen- 
dence and magnitude of the changes in the phonon fea- 
tures. The right panel presents the same comparison for 
the lone phonon. This figure shows that, even though 
Eq.[T] neglects anharmonic effects, its outcoming parame- 
ters are representative of the physical properties of MnF2. 
We also checked that comparison of the Lorentz fits to 
other methods such as four-parameter simulations^^ and 
multi-oscillator fit 9^ give the same results. Henceforth 
we will discuss our results using the Lorentz oscillator 
model alone. 

We performed ab-initio calculations to predict the val- 
ues of infrared active phonon frequencies and oscillator 
strengths. The results of these calculations for the in- 
frared active modes are shown as dashed lines in both 
panels of Fig. O For reference, the results for all modes 
are given in the appendix. The ab-initio frequencies 
are between 2 and 12% higher than the measured ones. 
Except for the very weak Eu2 phonon, the predicted 
strenghs are underestimated by 4% to 35%. Neverthe- 
less, the calculated spectra describe very well the overall 
phonon response. Table |I] summarizes the fitting param- 
eters produced by Eq. [T]at 5 K, 100 K and 300 K. It also 
shows the ab-initio results for Qto and As in the FM 
and AFM configurations defined in Sec. [Ill 

X-ray diffraction results are presented on Fig 2] where 
we plot the relative change of both lattice parameters 
as a function of temperature. We also plot in this fig- 
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FIG. 3: (color online) The left panel shows €2 spectra ob- 
tained by Kramers-Kronig transformation (solid curves) for 
phonon Eui and the right panel for phonon A2u- For clar- 
ity, each S2 curve was normalized by its maximum value then 
shifted vertically by a value proportional to its temperature 
(the curves are placed so that their maxima coincide with the 
temperature scale shown). The squares are ^to frequencies 
obtained from the Lorentz fits and are shown in a temperature 
vs. frequency plot. In a perfect agreement between Lorentz 
and Kramers-Kronig analysis, each square would sit exactly 
on the maximum of a curve. 
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FIG. 4: (color online) Relative changes in a and c lattice pa- 
rameters measured by x-ray diffraction. We also show the 
resulting relative variation of the unit-cell volume (open tri- 
angles) . 



ure the relative change in the unit cell volume given by 
AV/V = 2 (Aa/a) + (Ac/c) + [1 + (Ac/c)] {Aa/af + 
2 (Aa/a) (Ac/c). The c lattice parameter, as expected, 
increases monotonically with temperature. The AFM 
transition has a clear effect on this parameter which 
shows a kink at Tn- The most striking feature of this 
plot is the anomalous behavior of the a lattice parameter 
which increases with decreasing temperature and shows 
absolutely no feature at T^. Despite the anomalous a 
axis behavior, the total volume of the lattice still in- 
creases with increasing temperature. It is also worth 
mentioning that no change was observed on lattice pa- 
rameters when a 3 T magnetic field is applied. 
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TABLE I: Lorentz fit parameters at 5 K, 100 K and 300 K and ab-initio results for tfie pfionon frequencies and oscillator 
strengths. Units for Qto and 7 are cm~^. FM and AFM configurations are defined in Sec. [HI Fitted values for £00 are 2.16 
for the a6-plane and 2.25 along the c axis. 
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IV. DISCUSSION 

As a general rule, a crystalline lattice will contract 
upon cooling the sample. At a magnetic transition fur- 
ther modifications of the crystalline parameters are in- 
duced by the magnetic energy change. Indeed, following 
the classical work of Baltensperger and Helman^^ one 
can derive the modification of the lattice parameters in- 
duced by the magnetic ordering (at a given temperature 
T below Tn) by the minimization of the free energy (J^). 
Below Tat, the entropy term of the free energy is negli- 
gible and can be written as the sum of the magnetic 
(Em) and the elastic {Se) energies. A Taylor's expansion 
of as a function of the lattice parameters a and c yields 



J'(a, c) f (a, c) = Sm{cl, c) + fe(«, c) 



t + -—da + —dc - 
oa oc 



(2) 



- (aAr,C7v) 



where (a at, cat) are the lattice parameters at T < T/v in 
the antiferromagnetic phase. Naming ao,co the lattice 
parameters at the same temperature in an hypothetic 
paramagnetic phase, one gets 
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(4) 



The elastic energy per unit cell can be estimated follow- 
ing the classical way of Ref. [s^: 



_V(AV\' 



(5) 



where V is the unit cell volume and the compressibility 
coefficient. 



TABLE II: Local quantification axes of the two Mn atoms in 
the unit cell. The unit cell framework is expressed in terms 
of its main a, b and c directions, whereas x, y and z refer to 
the local Mn quantification axes. 



Local axes 

X 

y 

z 



Direction in the unit cell framework 
Mm (0,0,0) Mn2 (1/2,1/2,1/2) 



{-a + b)/V2 
c 



(a + b)/V2 
—c 

{-a + b)/V2 



Figure [5] and Tab. |TT1 show the definitions for the lo- 
cal quantification axes (x, y and z) of each Mn atom in 
the unit cell framework. The local Mn (x, y) plane is de- 
fined by the in-plane (equatorial) fluorine atoms of the 
distorted octahedra (shown as the thick line in Fig. [5|), 
which have the strongest metal-ligand interactions. Note 
that the two manganese atoms of the unit cell have dif- 
ferent local quantification axes. Using this definition, the 
magnetic energy can be expressed as 



(6) 



It is important to notice that in the case of MnF2, 
where the interactions are antiferromagnetic, both J and 
(5'Mni • ^S^Mns) ^rc negative. J can be decomposed into its 
main superexchange coupling. Using the local quantifica- 
tion axes defined in Fig.[5]and Tab.HIl we determined the 
superexchange coupling paths and their associated per- 
turbative contributions pictured in Tab. Hill The sum of 
these terms leads to: 

J = — exp [—2\{da + de)] X 

(Aq + A2 cos^a + cos^a + Aq cos^a) . (7) 

da and de represent the apical and equatorial Mn-F dis- 
tances, a is the angle Mni-F-Mn2. Although the ab- 
solute values of the Ai coefficients are difficult to deter- 
mine, we can estimate their relative values. We found 
A4 -Aq ^ 10^ A2 ^ 10^ Aq. Therefore the last two 
terms in Eq|71 dominate the value of J. 

We can use the above equations to obtain a rough eval- 
uation of the derivatives in Eq. |4j We find that the first 
derivatives are of the same sign and order of magnitude. 



TABLE III: Main AFM coupling paths between adjacent Mn 
atoms. The Mn d orbitals are written in terms of their local 
axes (see Tab. Ulland Fig. [5]) and the F p orbitals in terms of 
Mn2 local axes. The vertical and horizontal axes in pictures 
correspond to the d-\-b and —c lattice directions, respectively. 
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FIG. 5: (color online) Structure of MnF2 highlighting the two 
possible orientations for the MnFe coordination octahedra. 
The thick dark lines connect the in-plane F atoms. The local 
quantification axes orientations for Mni and Mn2 atoms are 
shown as black arrows. The lattice coordinates framework is 
defined by the arrows labeled a, b and c. 



This is also the case for d'^E/dadc and d'^E/dc^. On 
the contrary, d'^£/do? is nearly an order of magnitude 
larger than the other second derivative and of opposite 
sign. The reason is that only in the latter term magnetic 
and elastic contributions do not contribute with oppo- 
site signs. These estimations lead to a very weak value 
for Aa, while Ac is much larger and of opposite sign, in 
agreement with the x-ray measurements shown in Fig. |4j 

As the atoms get closer to each other, the interaction 
forces between them get stronger and, in the harmonic 
approximation, so does the spring constant {n) between 
ions. As phonon frequencies follow ~ \/~i^JJji {iJi being a 
reduced mass), in a first approximation we expect them 
to follow the thermal evolution of the unit cell volume. 
Figure m shows that the volume of the MnF2 lattice in- 
creases with temperature. Therefore, we expect phonon 
frequencies to decrease with increasing temperature. Fig- 
ures [6] and [7] show the thermal evolution of phonon fre- 
quencies for MnF2. Phonons Eu2 and Eu^ show the ex- 
pected frequency softening with increasing temperature 
over the whole measured range. Phonon also behaves 
conventionally, but only above Tn- 

Phonon Eui^ on the other hand, has an anomalous 
thermal evolution as its frequency softens with decreas- 
ing temperature. This behavior is analogous to that re- 
ported for the Big mode seen by Raman spectroscopy in 
MnF2 (Ref. 14) and other rutile type fluorides such as 
ZnF2,^^ FeF2,^° NiF2,^^ and MgF2.^^ This mode seems 
to be more influenced by the a parameter rather than the 
whole unit cell volume. The unconventional lattice pa- 
rameters thermal evolution and this anomalous phonon 
behavior are strong indications of large lattice instabili- 
ties. 
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FIG. 6: (color online) Evolution of the ab plane phonon fre- 
quencies with temperature. The red arrows indicate the ab- 
initio predictions for the change in frequencies (numerically 
indicated as the 6(jJ values) upon formation of an antiferro- 
magnetic order. The shaded area is the antiferromagnetic 
phase. For each vibration mode, we represent the atomic 
motions of manganese and fluorine ions. The dotted lines are 
guides to the eye showing the high temperature extrapolations 
of the phonons behavior in the absence of magnetic ordering. 



The existence of these lattice instabilities is empha- 
sized by Raman scattering under pressure. In ZnF2 
the softening of the Big Raman mode with decreasing 
temperature is also observed by applying pressure at 
room temperature. This pressure induced phonon soft- 
ening is a precursor of a phase transition from rutile 
to a CaCl2 structure at 4.5 GPa.^^ Brillouin scatter- 
ing measurements^^ showed that MnF2 also undergoes 
a phase transition at 1.49 GPa, whose nature was incon- 
clusive. Since the transition in ZnF2 is accompanied by 
an orthorhombic distortion of the lattice, we can inter- 
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FIG. 7: (color online) Evolution of the c axis phonon fre- 
quency with temperature. The red arrow indicates the value 
of the change in frequency when the antiferromagnetic order 
appears, as predicted by ab-initio calculations. The shaded 
area is the antiferromagnetic phase. The atomic motions of 
manganese and fluorine ions are given for this vibration mode. 
The dotted line is a guide to the eye showing the high tem- 
perature extrapolation of the phonon behavior in the absence 
of magnetic ordering. 



pret the Eui mode softening in MnF2 as a precursor of a 
(incipient) phase transition. Further insight on this tran- 
sition is given by measurements of the elastic properties 
in rutiles. Melcher^^ showed that the C44 elastic constant 
in MnF2 has one anomaly at Tn and another, unrelated 
to magnetism, at lower temperatures. Rimai^^ observed 
a similar low temperature accident in the non magnetic 
rutile ZnF2 and concluded that it was representative of 
a structural instability compatible with ferroelectricity. 
Hence, although only new measurements will be able to 
settle the issue, it is reasonable to expect that the high 
pressure transition in MnF2 is of (incipient) ferroelectric 
character. 

Figures [6] and [3 show that at the AFM transition, there 
is a clear renormalization of the phonon spectra. Mode 
Eu3 shows an additional increase in its frequency below 
Tat. This is an effect compatible with the magnetostric- 
tive kink observed in the lattice volume, which decreases 
faster below T^. There is no noticeable change in the be- 
havior of phonon Eu2' A puzzling behavior is observed 
for phonons Eui and A2U' The phonon softening instabil- 
ity observed in mode Eui stops at Tn and the frequency 
of this phonon increases notably with the appearance of 
the AF phase. Conversely, the thermal evolution of mode 
reverses, and this phonon shows a frequency soften- 
ing below Tat, which is at odds with the lattice volume 
evolution. As no lattice parameter show a sign change 
in their thermal behavior at Tat, the sign change in the 



8 



slope of the thermal evolution of these two phonon modes 
is likely due to effects other than conventional magne- 
tostriction, such as exchange between ions. 

To grasp further insight on this phonon behavior, we 
can look into ab-initio results. The paramagnetic phase 
is not accessible to ab-initio calculations. We consid- 
ered that the average of calculations in the AFM and 
FM configurations are representative of the paramagnetic 
phase. The difference in the phonon frequencies between 
the AFM and "paramagnetic" calculations are indicated 
by the vertical arrows in Figs. [6] and [3 It is remarkable 
that not only do ab-initio calculations predict the correct 
magnitude for the phonon frequency shifts, but also they 
find the proper sign for these changes. It is important 
to remark that ab-initio calculations are done at K. 
Hence, its predictions for phonon changes are not due 
to lattice modifications but are representative of a direct 
coupling of the phonons to the magnetic ordering, with- 
out the need of mediation by elastic lattice distortions. 

Ab-initio calculations also yield the normal mode 
eigenvectors. One can thus analyze the effect of the dis- 
placements associated with these phonon modes on the 
magnetic coupling constant. The infrared active normal 
modes are shown in Figs. [6] and [71 The thermal evolu- 
tion of all phonons in the paramagnetic phase, including 
the anomalous softening when lowering the temperature 
of the Eui mode is probably dominated by the elastic 
energy. Indeed, the same behavior is observed in the 
Raman spectra of the Big mode of MgF2 (Ref. 1351) and 
the infrared response of the A2u mode of ZnFs (Ref.liO), 
both non magnetic fluorides, isostructural to MnF2. 

In order to understand the sign change in the thermal 
evolution of phonons Eui and below Tat, we have to 
look at the free energy = £ — TS. In the AFM phase, 
the dominant term in the changes of the free energy of the 
system is the magnetic energy £m, given by Eq. [6l Thus, 
minimization of goes through the minimization of Sm- 
Let us look into the contribution coming from phonons 
Eui and A2u to this term. The Eui mode corresponds to 
a large displacement in the a -\- b and a — b directions of 
the F atoms bridging two MnFe octahedra. Let us define 
as the magnetic coupling between the Mn2 atom and 
the Mni atoms located at the crystallographic positions 
(0,0,0) and (0,0,1); J^_^^ as the coupling between Mn2 
and the Mni atoms located at (1, 1,0) and (1, 1, 1); and 
as the coupling between Mn2 and the Mni atoms lo- 
cated at (1, 0, 0), (1, 0, 1), (0, 1, 0) and (0, 1, 1). Under a 
Eui phonon displacement of amplitude the magnetic 
energy is modified as follows 



8 J(0) ■ 



u=0 



{Suni ■ S'Mns)- (§) 



J(0) is the value of the magnetic exchange with the atoms 
in their equilibrium position. We should stress, again, 
that both J(0) and the spin correlation terms are nega- 



tive because of the AFM interaction. Using the pertur- 
bative expressions of Tab. Illll one can show that does 
not depend on the Eui displacement that dJ^/du = 
-dJ^j^^^/dUi and that d^J^jdu^ = d^J^^-^jdu^ > 0. In 
the AFM phase the magnetic energy must be minimized, 
and thus, the modulus of the first factor of Eq. [8] must 
increase. As the second derivative and J(0) have oppo- 
site signs, increasing the modulus of the first factor im- 
plies that the oscillation amplitude u must decrease. In a 
equivalent view, one can note that the second derivative 
term acts as a correction to the full harmonic potential 
of the phonon. For phonon Eui it introduces a positive 
correction to the effective phonon spring constant, i.e. 
the phonon hardens, as observed in Fig. [6l 

A similar analysis also explains why the A2u mode soft- 
ens when the temperature decreases in the AFM phase. 
The A2u mode corresponds to displacements with op- 
posite signs of the Mn and F atoms along the c direc- 
tion (see Fig [7]). We define as the magnetic coupling 
between the Mn2 atom and the Mni atoms located at 
(0, 0, 0), (1, 1, 0), (1,0, 1), (0, 1, 1); and as the coupling 
between the Mn2 atom and the Mni atoms located at 
(0, 0, 1), (1, 1, 1), (1, 0, 0), (0, 1, 0). Under a A2u displace- 
ment of amplitude the magnetic energy is modified 
as 



£m 



[4 J5(v) + 4 Jc'(v)] (^Mni • ^Mn; 

8J(0) ' 



{Syini • SMn2)' (9) 



Using again a rough evaluation of the perturbative ex- 
pression of J one finds that dJ^/dv = —dJ^/dv^ and 
d^J^jdv^ = d'^J^/dv'^ < 0. Because the second deriva- 
tive and J(0) here have the same sign, an increase of the 
displacement amplitude v diminishes the value of £m- In 
the full harmonic potential perspective, here the correc- 
tion to the effective phonon spring constant is negative, 
leading to the phonon softening, in agreement with Fig. [71 
Our infrared data allow us to link the magnetic or- 
dering to changes in the dielectric constant of MnF2. 
The optical dielectric constant is given by £:(0) = £00 + 
Asfc. In the absence of any microwave excitations, it 
should be equal to the static value obtained from elec- 
trical measurements. Our reflectivity data do not show 
any extra excitation below the phonons down to 10 cm~^ 
(300 GHz). Figure [HI compares the static dielectric con- 
stants, measured by Seehra^^ at 10 kHz, with the zero 
frequency limit of our optical data. The left panel shows 
the a6-plane dielectric constant, and the right panel the 
c axis response. The sohd squares are the "static" val- 
ues from Seehra^^ (vertically shifted by -1, for clarity) 
and the open circles are the optical dielectric constants. 
On the left panel, we also show the thermal variation of 
the oscillator strength for the lowest energy Eui mode 
(solid triangles, right-hand scale). Note that both scales 
in this figure cover the same variation in values of e. It 
is then clear that all the temperature variations in the 
optical dielectric constant come from the Eui phonon. 
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FIG. 8: (color online) The left panel shows the ab plane di- 
electric constant £(0) measured utilizing dielectric techniques 
by Seehra^^ (open squares) and zero frequency optical ex- 
trapolation (solid circles). The right panel shows the same 
quantities for the c axis. For clarity, the dielectric measure- 
ment was shifted by -1. The blue arrows are the ab-initio 
predictions for changes Se in ^(0) upon appearance of the an- 
tiferromagnetic order. The right-hand scale in the left panel 
shows the Eui contribution to the dielectric constant (solid 
stars) and the red arrow is the ab-initio prediction for the 
change in A£{Eui). Note that right and left scales in this 
panel span the same range in values of e, indicating that all 
the in-plane changes come from the Eui mode. 



This figure also shows that, within 15%, the static dielec- 
tric constant is build from phonon and higher frequency 
excitations. The static and optical dielectric constants 
agree very well below Tn. In the paramagnetic phase, 
there is a growing difference between these two values. 
Along the c-axis both values increase with temperature 
but the static dielectric constant increases faster. In the 
a6-plane the difference is also qualitative. The static di- 
electric constant increases with temperature, whereas the 
optical value decreases. These discrepancies may be due 
to a difference in the quality of the sample, but are more 
likely linked to ionic conduction contributions to the di- 
electric constant in the radiofrequency range, a common 
effect on fluorides. ^^^^^ Nevertheless, in both polariza- 
tions, £(0), measured by both techniques, shows a kink 
at T/v- 

We can estimate the ab-initio zero frequency dielectric 
constant from the calculated oscillator strengths. We 
made the same comparison between a "paramagnetic" 
and AFM configurations, as we did for the phonon fre- 
quencies. The predicted changes are shown as the verti- 
cal arrows in Fig. [8l The predicted ab-initio jump upon 
magnetic ordering is of the same order of magnitude as 
the values measured. More important, we find the cor- 
rect sign for the changes in both cases. The AFM order 
tends to decrease the dielectric constant in the a6-plane 
while it is increased along the c-axis. 

The panorama in MnF2 points towards instabilities 
that favor a ferroelectric order. Ferroelectricity is not 
a common property of fluorides. However, several com- 
pounds (ZnF2, NiF2 and FeF2) show lattice instabilities 
that indicate the appearance of an incipient ferroelec- 
tric order. Our x-ray data indicate that lattice insta- 
bilities are indeed present in MnF2. Infrared measure- 



ments show that the c-axis phonon softens upon forma- 
tion of the antiferromagnetic order. This softening goes 
together with an increase of the dielectric constant along 
this direction. These effects are indicative of a possi- 
ble incipient ferroelectric behavior in MnF2. In this sce- 
nario the differences between static and phononic dielec- 
tric constants observed in Fig. [8] above 100 K could also 
be related to zero frequency ferroelectric fluctuations. 
However, dielectric constant measurements at high fre- 
quencies as well as phonon dynamic measurements un- 
der pressure and/or magnetic field are necessary to con- 
firm (or refuse) this picture. In any case, understanding 
the phonon changes induced by the magnetic ordering in 
MnF2 should bring new insights on the magneto-electric 
coupling in multiferroic materials. 



V. CONCLUSIONS 

In this paper we presented a detailed temperature de- 
pendent infrared study of the phonon renormalization ob- 
served at Tn in MnF2 . We showed that the phonon tem- 
perature dependence and the lattice parameter changes 
across the antiferromagnetic transition are well repro- 
duced by ab-initio calculations, implying that the mag- 
netic order dominates the changes observed at Tat. We 
find that phonons along the D4 axis and the tetragonal 
ab plane have opposite changes at Tn as predicted by 
the first principles results. Our results show that the di- 
electric constant of MnF2 is mostly from phonon origin. 
Relaxation effects on the ab plane contribute strongly to 
the dielectric constant in the paramagnetic phase. The 
large change in the lattice parameters with temperature 
and the phonon softening in the antiferromagnetic phase 
suggest that MnF2 distortions induced by the magnetic 
order are compatible with the ferroelectric instabilities 
observed in Ti02 and FeF2 and other fluorides. We 
suggest that the phase transition observed at 1.49 GPa 
(Ref. [37I ) could be of ferroelectric order and so, MnF2 
would qualify as an incipient multiferroic. 
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Appendix 

In this paper, we limited our discussion of ab-initio 
calculations to the four infrared active modes only. For 
completeness, Table IIVI gives the results of calculations 
for all modes (infrared, Raman and silent) in both FM 
and pseudo-AFM configurations. These calculations are 
compared to the experimental values determined at 5 K 
for the infrared (this work) and Raman active modes 
(Ref.[ll. 



TABLE IV: Ab-initio F point optical phonons frequencies (in 
cm~^) for the FM and pseudo-AFM phases. The experimen- 
tal values were taken at 5 K from the present work (IR) and 
from Ref. 14 (Raman) 



ab-initio Experimental 
Symmetry Activity FM AFM This work Ref. JA 



Big 


Raman 


96.2 


97.8 




56.1 






136.2 


136.4 






Eu 


IR 


175.2 


178.2 


158.4 




A2g 




238.9 


241.1 






Eg 


Raman 


257.8 


257.2 




246.7 


Eu 


IR 


265.5 


265.4 


259.8 




A2u 


IR 


315.1 


311.0 


292.5 




Blu 




346.7 


344.4 






Alg 


Raman 


365.8 


368.7 




346.8 


Eu 


IR 


387.7 


394.6 


367.0 




B2g 


Raman 


490.6 


491.6 




480.5 
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